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All possible linearly independent local conservation laws for n-dimensional difFusion-convection 

f— ^ I equations Ut — {A(u))ii + {B^{u))i were constructed using the direct method and the composite 

f^ ■ variational principle. Application of the method of classification of conservation laws with 

CN ' respect to the group of point transformations [R.O. Popovych, N.M. Ivanova, J. Math. Phys., 

{3JQ. 2005, V.46, 043502 | |math-ph/0407008 1] allows us to formulate the result in exphcit closed form. 
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^ ' Action of the symmetry groups on the conservation laws of diffusion equations is investigated 

•^ ■ and generating sets of conservation laws are constructed. 

CN ■ 1 Introduction 

> 

L~ ■ One of the most important applications of group analysis is the construction of conservation 

(^ , laws of (systems of) differential equations. They provide information on the basic properties of 

^T I solutions of differential equations. The famous laws of conservation of energy, linear momentum 

\^ ' and angular momentum are important tools for solving many problems arising in mathematical 

^^ . physics. Knowledge of conservation laws is important for numerical integration of PDEs. Inves- 

tigation of conservation laws of Korteweg-de Vries equation became a starting point of discovery 
Oh' of new approaches to integration of PDEs (such as Miura transformations, Lax pairs, inverse 

r~| ■ scattering, bi-Hamiltonian structures etc.) Existence of the 'sufficient number' of conservation 

j^ I laws of (systems of) PDEs is a reliable indicator of their possible integrability 

In view of the generalized Noether theorem, there exists a one-to-one correspondence between 
the non-trivial generalized variational symmetries of some variational functional and the non- 
trivial conservation laws of the associated Euler-Lagrange equations, and any such symmetry 
/\ ' is a generalized symmetry of the Euler-Lagrange equations. Thus, e.g., conservation of energy 

j^ ■ arises from the invariance of the corresponding variational problem with respect to the group 

of time translation. The Noether approach reduces construction of conservation laws to finding 
symmetries for which there exist a number of well-developed methods. However, this approach 
can be applied only to Euler-Lagrange equations that form normal systems and admit symmetry 
groups satisfying an additional "variational" property of leaving the variational integral invariant 
in some sense [19]. The latter requirements lead to restriction of class of systems that could be 
investigated in such way. 

At the same time, the definition of conservation laws itself yields a method of finding con- 
servation laws. This method is called direct, and it is possible to distinguish four its versions, 
depending on the way of taking into account systems under investigation. (See e.g. [2,3,19,20,24] 
for description and detailed comparison of different methods of finding conservation laws.) In 
the present paper we use its most direct version based on the immediate solving of determining 
equations for conserved vectors on the solution manifolds of investigated systems. 



Let us note that there exist other approaches for construction of conservation laws which 
differ from the Noether or direct ones, are based on expfoitation of symmetry properties of 
differential equations and can be applied to non-Lagrangian systems. Thus, W.I. Fushchych 
and A.G. Nikitin [10] proposed to directly calculate bilinear combinations of solutions of motion 
equations, which are conserved in time by virtue of symmetries of these equations. 

Another useful approach for construction of conservation laws of non-Lagrangian systems is 
so-called composite variational principle [4,23]. Composite principle means that in addition to 
the original variables of a given system, one should introduce a set of adjoint variables in order 
to obtain a system of Euler-Lagrange equations for some variational functional. To the best of 
our knowledge such method was firstly developed and applyed for finding conservation laws of 
the linear one-dimensional heat equation [4]. Recently this approach was rediscovered and used 
by N.H. Ibragimov and T. Kolsrud for investigation of conservation laws of one-dimensional 
diffusion equations [14]. 

The aim of this work is to construct all possible linearly independent local conservation laws 
of n-dimensional diffusion-convection equations of the form 

ut = {A{u)h + {B\u))i, (1) 

where A{u) and B^{u) are arbitrary smooth functions, A^ ^ 0, u = u{x) = u{x^, . . . ,x^), 
Ui = du/dx^, Uij = d'^u/dx'^dxK Here and below we use the summation convention for repeated 
indices, i = l,...,n. This equation, often called the Richard's equation, arises naturally in 
certain physical applications. Thus, for example, superdiffusivities of this type have been pro- 
posed [11] as a model for long-range Van der Waals interactions in thin films spreading on solid 
surfaces. This equation also appears in the study of cellular automata and interacting particle 
systems with self-organized criticality (see [6] and references therein). It describes a model of 
water flow in unsaturated soil [21]. 

The main tool of our investigation is the notion of equivalence of conservation laws with 
respect to equivalence groups which was introduced in [20]. Investigation of such equivalences 
allows us to simplify essentially technical calculations and to obtain results in closed explicit form. 
Note that for wide classes of nonlinear equations, this is the only way to obtain the complete 
description of spaces of conservation laws (see, e.g., [16,20] for detailed analysis of application of 
such equivalences for finding conservation laws of different classes of one-dimensional diffusion- 
convection equations). 

Conservation laws of equations (^ were studied firstly for its one-dimensional case. Thus, 
A.H. Kara and F.M. Mahomed [17] constructed the first-order local conservation laws of equa- 
tions (^ with n = 1 and found the bases of such conservation laws with respect to the corre- 
sponding symmetry groups. In [15,20] the complete classification of local and potential conser- 
vation laws for one-dimensional equations (^ with respect to the equivalence group is presented. 
V.A. Dorodnitsyn and S.R. Svirshchevskii [9] (see also [13, Chapter 10]) completely investigated 
the local conservation laws for one-dimensional reaction-diffusion equations ut = {A{u))xx+C{u) 
having non-empty intersection with class (^. Conservation laws of multidimensional reaction- 
diffusion equations ut = A{A{u)) + f{u) were found by Y.R. Romanovsky [22], who completely 
studied conservation laws for case A{u) = u™", and R.M. Cherniha [7], who considered the cases 
of power and exponential nonlinearities. 

In the present work we briefly review the necessary theoretical background including the 
notions of equivalence of conservation laws with respect to groups of equivalence and symmetry 
transformations (Section 2) introduced in our recent work [20], construct the complete equiv- 
alence group and exhaustively classify, with respect to this group, the local conservation laws 
of equations (^ (Section 3). Then, in Section 4, we apply the composite variational principle 



to the subclass of diffusion equations with zeroth convectivities, investigate the classical Lie 
symmetries of the obtained system of Euler-Lagrange equations, and use them to construct 
conservation laws of the systems by means of Noether approach. In such a way we found a set 
of local and nonlocal conservation laws for the diffusion equations, where the non-local variable 
is defined by the associated equation for the adjoint symmetry. Action of symmetry groups on 
the conservation laws of diffusion equations is investigated and generating sets of conservation 
laws are found in Section 4. 

2 Theoretical Background 

Let £ be a system L(x, ti(p)) = of / differential equations L^ = 0, . . . , L' = for m unknown 
functions u = (n^, . . . , u^) of n independent variables x = (x^, . . . , x"). Here U(p) denotes the set 
of all the derivatives of the functions u with respect to x of order not greater than p, including 
u as the derivatives of the zero order. Let £(fc) denote the set of all algebraically independent 
differential consequences that have, as differential equations, orders not greater than k. We 
identify £(;-) with the manifold determined by >C(fc) in the jet space J^''' . 

Definition 1. A conserved vector of the system C is an n-tuple F = {F^(x, u/.^\), . . . , F^{x, u/.f.\)) 
for which the divergence DivF := DiF^ vanishes for all solutions of C (i.e., DivFL= 0). 

In definition n and below Di = D,j.i denotes the operator of total differentiation with respect 
to the variable x*. The notation FL means that values of V are considered only on solutions of 
the system C 

Definition 2. A conserved vector F is called trivial if F* = F^ + F^, i = l,n, where F^ and 
F* are, likewise F*, functions of x and derivatives of u (i.e., differential functions), -F* vanish 
on the solutions of C and the n-tuple F = (F^, . . . , F") is a null divergence (i.e., its divergence 
vanishes identically). 

Definition 3. Two conserved vectors F and F' are called equivalent ii the vector- function F' — F 
is a trivial conserved vector. 

For any system C of differential equations the set CV(£) of conserved vectors is a linear 
space, and the subset CVo(£) of trivial conserved vectors is a linear subspace in CV(£). The 
factor space CL{C) = CY{C)/ CVo{^) coincides with the set of equivalence classes of CV(£) 
with respect to the equivalence relation adduced in definition |21 

Definition 4. The elements of Ch{C) are called conservation laws of the system £., and the 
whole factor space CL(£) is called the space of conservation laws of C. 

That is why we assume description of the set of conservation laws as finding CL{C) that 
is equivalent to construction of either a basis if dimCL(/^) < cxo or a system of generatrices 
in the infinite dimensional case. We will additionally identify elements from CL(i2) with their 
representatives in CV(£). Namely, we will understand a conservation law as a divergence ex- 
pression F'jF* = vanishing identically for all solutions of C. In contrast to the order rp of a 
conserved vector F as the maximal order of derivatives explicitly appearing in F, the order of a 
conservation law as an element J- from CL(/3) is called in.m{rF \ F € CV(>C) corresponds to J-}. 
Under the linear dependence of conservation laws, we understand the linear dependence of them 
as elements of CL{C). 

Definition 5. Conservation laws of a system C are called linearly dependent if there exists their 
linear combination having a trivial conserved vector. 



Let the system L be totally nondegenerate [19]. Then application of the Hadamard lemma 
to the definition of conservation law and integrating by parts imply that up to the equivalence 
relation of conserved vectors the left-hand side of any conservation law of L can be always 
represented as a linear combination of the left-hand sides of independent equations from L with 
coefficients A'^ being functions on a suitable jet space J^^': 

BwF = X''L''. (2) 

Here the order k is determined by £. and the allowable order of conservation laws, /i = 1,1. 

Definition 6. Formula ^ and the /-tuple A = (A^, . . . , A ) are called the characteristic form 
and the characteristic of the conservation law Div F = correspondingly. 

The characteristic A is trivial if it vanishes for all solutions of C. Since C is nondegenerate, 
the characteristics A and A satisfy © for the same F and, therefore, are called equivalent iff 
A — A is a trivial characteristic. Similarly to conserved vectors, the set Ch(>C) of characteristics 
corresponding to conservation laws of the system £ is a linear space, and the subset Cho(£) of 
trivial characteristics is a linear subspace in Ch(£). The factor space Chf(£) = Ch(£)/ Cho(i2) 
coincides with the set of equivalence classes of Ch(£) with respect to the above characteristic 
equivalence relation. 

Using properties of total divergences, we can eliminate the conserved vector F from @ and 
obtain a condition for the characteristic A only. Namely, a differential function / is a total 
divergence, i.e. / = Div F for some n-tuple F of differential functions iff E(/) = 0. Hereafter 
the Euler operator E = (E^, . . . , E™) is the m-tuple of differential operators 



E« = (-Z))-5„., a = l,m, 

where a = (ai, . . . ,a„) runs the multi-indices set (ajGNU{0}), (— -D)" = (— Di)"i . . . {—D, 
Therefore, action of the Euler operator on ^ results to the equation 



mj 



E{X^Ln = Dl{L) + DUX) = 0, (3) 

which is a necessary and sufficient condition on characteristics of conservation laws for the sys- 
tem C The matrix differential operators D^ and D2 ^^^ the adjoints of the Frechet derivatives D;^ 
and D^, i-e.. 

Since D^(^) = automatically on solutions of C then equation ^ implies a necessary condition 
for A to belong to Ch{C): 

Di(A)|^=0. (4) 

Condition ^ can be considered as adjoint to the criteria D/^(?/)L= for infinitesimal invariance 
of C with respect to evolutionary vector field having the characteristic r/ = (r/^, . . . ,??*"). That 
is why solutions of Q are sometimes called as cosymmetries [5] or adjoint symmetries [3]. 

Let now C he a system of Euler-Lagrange equations E L = for the variational functional 
^ = f L{x,U(^p))dx. 

Definition 7. Generalized vector field Q = S,^di + r]°'du'^ is called variational symmetry of 
functional ^ = J \_{x, U(^p-^)dx iff there exists a tuple of m differential functions {B^, . . . , B"^), 
such that (5(p)(L) + L-DiC* = DiB'^ for all x and u. 



It can be proved [19] that any variational symmetry of the variational problem is a symmetry 
of the corresponding Euler-Lagrange equations. (In general, the inverse statement is not true.) 
Thus, to find all variational symmetries of the given variational functional, it is enough to con- 
struct symmetries of the corresponding Euler-Lagrange equations and then to check additional 
criterion of variational invariance given in definition d 

Variational symmetry Q is trivial if its characteristic Q[u\ = [r]^ — ^^uj, . . . ,7]"^ — C'^^T) 
vanishes on the solution manifold of the Euler-Lagrange equations. Two variational symmetries 
Q^ and Q^ are called equivalent if Q^ — Q^ is a trivial symmetry. 

There exists one-to-one correspondence between the non-trivial generalized variational sym- 
metries of the variational functional and the non-trivial conservation laws of the associated 
Euler-Lagrange equations, and this correspondence is established by the generalized Noether 
theorem [19]: 

Theorem 1. Generalized vector field Q determines the group of variational symmetries of func- 
tional ^ = J \_dx iff its characteristic is a characteristic of a conservation law of the associated 
Euler-Lagrange equations E(L) = 0. In particular, if ^ is a nondegenerate variational problem, 
there exists a one-to-one correspondence between the equivalence classes of nontrivial conserva- 
tion laws of the Euler-Lagrange equations and equivalence classes of variational symmetries of 
the functional. 

In most of cases to describe completely the space of conservation laws of (a class of) systems 
of differential equations, one should take into account symmetry transformations of a system or 
equivalence transformations of a whole class of systems [20] . 

Proposition 1. Any point transformation g maps a class of equations in the conserved form 
into itself. More exactly, the transformation g: x = Xg{x,u), u = Ug{x,u) prolonged to the 
jet space J^^' transforms the equation DiF^ = to the equation DiF^ = 0. The transformed 
conserved vector Fg is determined by the formula 

Dx Xi ■ 1 

Fl{x,ui^)) = ' F\x,U(r.)), i.e. Fg{x,u,^)) = ——-{Dxx)F{x,u,^)) (5) 

in the matrix notions. Here |-Da;5;| is the determinant of the matrix D^x = {Dx Xi). 

Note 1. In the case of one dependent variable {m = \) g can be a contact transformation: 

Jb Jb ri\ Jb % Ub. 



g\ 



i(i)), 11(1) = 'Wg{i) (2^)^(1))- Similar notes are also true for the statements given below. 

Definition 8. Let G be a symmetry group of the system C. Two conservation laws with the 
conserved vectors F and F' are called G -equivalent if there exists a transformation g & G such 
that the conserved vectors Fg and F' are equivalent in the sense of definition 13 

Proposition 2. If system £. admits a one-parameter group of transformations then the infinites- 
imal generator Q = ^*5j + ri°'du<^ of this group can be used for construction of new conservation 
laws from known ones. Namely, differentiating equation with respect to the parameter e and 
taking the value e = 0, we obtain the new conserved vector 

F^ = -Q^r)F' + [D,e)F^ - iDje)F\ (6) 

Here Q(^r) denotes the r-th prolongation of the operator Q. 

Therefore, for any system of differential equations we can introduce a generating set of conser- 
vation laws in such a way that all inequivalent conservation laws can be obtained from this set by 
multiple using the symmetry transformations and taking linear combinations. Such generating 
sets, called also L-bases, were firstly described in [18] (see also [12]). 



Lemma 1. Any point transformation g between systems L and L induces a linear one-to-one 
mapping g^, between the corresponding linear spaces of conservation laws. 

Consider the class C\s of systems L(x, n(p), 9{x, U(^p-^)) = parameterized with the parameter- 
functions 9 = 6(x,U(^p)). Here L is a tuple of fixed functions of x, Uf^p-^ and 9. 6 denotes the 
tuple of arbitrary (parametric) functions 9{x,U(^p^) = (0^(x,n(p)), . . . ,9^{x,U(^p-j)) satisfying the 
additional condition 5(j;, U(p), 0(q)(x,U(p))) = 0. 

Let P = P{L, S) denote the set of pairs each from which consists of a system from C\s and a 
conservation law of this system. Action of transformations from an equivalence group G^ of the 
class C\s together with the pure equivalence relation of conserved vectors naturally generates 
an equivalence relation on P. Classification of conservation laws with respect to C will be 
understood as classification in P with respect to the above equivalence relation. This problem 
can be investigated in the way that it is similar to group classification in classes of systems of 
differential equations. Specifically, we firstly construct the conservation laws that are defined 
for all values of the arbitrary elements. Then we classify, with respect to the equivalence group, 
arbitrary elements for each of the systems that admits additional conservation laws. 

For more detail and rigorous proof of the correctness of the above notions and statements 
see [20]. 

3 Classification of conservation laws 

To classify the conservation laws of equations from class (^ we have to start our investigation 
from finding its equivalence transformations. Application of both the direct and infinitesimal 
methods to class (^ allows us to construct the complete equivalence group C. The following 
statement is true. 

Theorem 2. Any transformation from C^quiv j^^^ ^^g form 
i = e^t -\- El, X* = e^fiijx^ + ent + e2i, u = equ + £3, 

A = E^'^eIA-\-E8, B' = E^'^E5HijB^ -E7iU-\-Egi, (7) 

where e's and fi's are constants, e^e^Eq y^ 0, M = {jiij) € SO{n), i.e. M is an arbitrary special 
orthogonal matrix: detM = 1, M = M~^ . 

Any conservation law of equation from class (^ has the form 

DtT{t, X, U(,)) + DiX\t, X, n(,)) = 0, (8) 

where T and X^ are often called conserved density and conserved fluxes correspondingly. 
Since (^ is an evolution equation of the second order, then up to the usual equivalence re- 
lation of conservation laws, any conservation law of (^ is equivalent to one with the conserved 
density T = T{t, x, u) and the fluxes X* = X*(t, x,u,ui, . . . , Un) [12]. 

We search (local) conservation laws of equations from class (^ , applying the modification of 
the direct method, which was proposed in [20] and is based on using the notion of equivalence 
of conservation laws with respect to a transformation group. First, expand total differentiation 
operators in ((H)) on the solution manifold of (^: 

Tt + Tu{AuUii + AuuUiUi + Blm) + Xl + Xlm + X^Mij = 0, 

and split the obtained expression with respect to the unconstrained variables Uij. Coefficients 
oiuii and Uij, i / j give 

X' = -TuAuUi + L^l'^\t,x,u)uj + R\t,x,u), 



where U^'^'{t,x,u) = —U^''^\t,x,u). Considering the conserved vector with the fluxes X* = 
X* — DjU^'^^ = —TuAuUi + R!'{t,x,u), which is equivalent to the initial one, without loss of 
generality we can assume that L^*'-'! = 0. Splitting the rest of the conservation law with respect 
to different powers of Uj yields 

Tun = 0, -T.uAu + TuBl, + Ri = 0, Tt + B^ = 0. 

It follows that up to the usual equivalence relation of conserved vectors, the conserved density 
and fluxes in (jH)) have the form 

T = a{t,x)u, X^ = —aAuUi + OiA — aB^ , (9) 

where A, B^ and a(t,x) satisfy the classifying equation 

atu + AaA - UiB^ = 0. (10) 

Let us note that function a{t,x) is a characteristics of conservation law (|HJ, and differentiat- 
ing (|l()j) with respect to u gives exactly an equation for the adjoint symmetry for equation (^. 

Solving ((TU)) up to equivalence transformations (O we obtain the complete classification of 
local conservation laws of diffusion-convection equation (^ . 

Theorem 3. A complete list of G^ -inequivalent equations (^ having nontrivial conservation 
laws of form Q with density and fluxes 

T = a{t, x)u, X* = —aAuUi + aiA — aB^, 

is exhausted by the following: 

1. \/A,B': a = l. 

2. yA, B^ = ... = B'' = 0, B''^^ = A, 

B '^ , • • • , -B" ^ (1, u, A) are linearly independent, < k < n : 
a = a(x\ . . . , x''+^), an + . . . + Okk + ak+i,k+i + afc+i = 0. 

3. yA, B^ = . . . = B'' = 0, B''^^, ...,B'^^{1,u,A) are linearly independent, 
< k <n: a = a(x\ ... ,x^), Aq = 0. 

4. A = u,B^ = ... = B^ = 0, B^^^, ... ,B'^ ^ {I, u) are linearly independent, 
< k < n : a = a{t, xi, . . . , Xk), at + Aa = 0. 

Together with values A and B^ we also adduce the constrains for characteristic a{t,x). 

Proof. Let us briefly sketch the proof of the theorem. There exist two essentially inequivalent 
cases of integrating equation (fTU]) : A E (1, u) and A (1, n). 

\. A ^ {l,u). Then, case B^ {l,u,A), 1 < i < n leads to the case 1 from the theorem 
list. Suppose now that there exist such j that B^ = cfA + cju + c'^ for all 1 < i < j. Then 
applying Galilean transformations x* = x* + cjt, B^ = B^ — c\u and translations of B^ from the 
equivalence group, without loss of generality we set c\ = c^ = 0. Thus, S* = c^A mod G~. If 
all c? = then we have B^ = 0, 1 < i < j and classifying equation splits (|in|) into the system 
at = ai = 0, i > j + 1, an = 0. Thus, we have case 3. Suppose now, that there exists such jo 
that C2° 7^ 0. Then up to the rotation transformations from the G~ we can set C2 = 0, for all 
l^^^j — 1)C-' = 1- Taking j = k + 1 and decomposing the classifying condition with respect 
to the linearly independent functions of u we obtain precisely the case 2. 

2. A e {l,u). Then A = u mod G~. If B' (1,^) for all i, then ai = 0, at = that yields 
a subcase of case 1. If there exist B^ E {l,u) then, applying the Galilean transformations and 
translations of B^ we deduce that B^ = mod C" that leads to case 4. D 



In the case of one space dimension n = 1 equations for characteristic a can be integrated 
exphcitly and we obtain simpler form of classification result [20] . 

Corollary 1. Any equation from class ^ with n = 1 has a conservation law of form Q where 

1. T = u, X = -AuUr,-B. (11) 

A complete list of G^ -inequivalent one- dimensional equations (jH)) having additional (i.e. linear 
independent with (|llj) ) conservation laws is exhausted by the following: 

2. VA, B = 0: T = xu, X = A- xAuU^, (12) 

3. yA, B = A: T = e^u, X = -e'^AuU^, (13) 

4. A = u, B = : T = au, X = axU — aux, (14) 

where x = x^ , X = X^ , B = B^ , a = a{t, x) is an arbitrary solution of the backward linear heat 
equation at + axx = 0. (Together with values A and B we also adduce complete lists of densities 
and fluxes of additional conservation laws.) 

Analysis of the classification given in theorem |31 and its corollary results in the following 
conclusions. The space of conservation laws for equations from class (^ in general case is one- 
dimensional and is generated by one with the conserved density T = u and fluxes X^ = —A^Ui — 
B^. Imposing restrictions on values of arbitrary functions A{u) and B^{u) leads to subclasses 
of Q with two-dimensional spaces of (local) conservation laws for the nonlinear equation in one 
space dimension and with infinite-dimensional spaces of conservation laws for linear equations 
and in the case n > 2. In this cases the space of conservation laws is parameterized with an 
arbitrary solution of a corresponding adjoint linear equation. 

4 Composite variational principle for finding conservation laws 

Another possible way of finding conservation laws of equations (^ is extension of the Lagrangian 
approach via introducing an 'associate' equation in order to obtain a system of Euler-Lagrange 
equation for some variational functional. An example for such an associate is the adjoint sym- 
metry equation [4,23]. Here we illustrate this approach on example of subclass 

ut = auii + auuf, (15) 

of equations (^ with zeroth convection terms B^ = 0. Here and below a = a{u) = Au{u). 

To find conservation laws of (|15j) from the composite variational principle we consider system 

ut = auuf + auii, vt + avii = (16) 

obtained by combining equation H15|) and the adjoint equation which is derivative of ()1U|) with 
respect to u. This is a system of Euler-Lagrange equations for the variational functional 

■^ = i-{utv-uvt)+auiVi]dx. (17) 

with Lagrangian L = ^{utv — uvt) + auiVi. Indeed, 

E" L = --wt + auViUi - -vt - Di{avi) = -{vt + avu), 

E" L = -Ut + -ut + Di{aui) = ut - {am + auu). 



Thus, we can use the Noether theorem for finding conservation laws of system ([T| 

In this work we restrict ourselves to the investigation of conservation laws corresponding 

to the Lie point symmetries of the Euler-Lagrange equations. (The complete classification of 

conservation laws of (|1H|) requires finding of all its higher order symmetries.) 

The Lie symmetries of system H15() will be classified up to the group G^ of equivalence 

transformations 

i=e4,t + ei, x^ = e^fiijX^ + e2i, u = eQU + e^, 

v = ev + a{x), a = e^ s^a, (18) 

where e, e's and /i's are constants, eiS^e^e ^ 0, M = (nij) S SO{n), i.e. M is an arbitrary 
special orthogonal matrix: detM = 1, M^^ = M~^, a{x) is arbitrary solution of the Laplace 
equation an = 0. This group is a natural extension of transformations (restricted to the case 
B^ = 0) onto the variable v. 

The classical Lie infinitesimal criterion yields the following system 

Ti = Tu = Ty = 0, ^i = C = 0, Vv = 0, C« = 0, 



a 



Vt - ar]jj = 0, Ci + aCjj = 0, 



aiCjj - 2mu) - 2aur]i - Ct = 0, aiCjj - 2C») + Ct = 0- 

of the determining equations for the coefficients of the infinitesimal symmetry generator 

Q = T{t,x,u,v)dt + C*(t,a;, u, v)di + ri{t, x,u,v)du + Cit,x,u,v)dy. 

Solving this system up to equivalence transformations (|18() we obtain the complete group clas- 
sification of system (|16|) . 

Theorem 4. The Lie algebra of the kernel of Lie invariance group of system ()16|) is the algebra 

^ker ^ ^g^^ q,^ x'dj-X^Oi, a{x)dy, vd^, 2tdt + x'di), 

where a{x) is arbitrary solution of the Laplace equation an = 0. 

A complete set of C^ -inequivalent systems (|16|) with the maximal Lie invariance algebras 
^max _^ ^kcr j^g exhausted by the ones with following parameter-functions a{u): 

1. a = e".- ^"^^^ = A^""' + {tdt - du), 

2. a = u^', (fi^ 0): A"^^^ = A^""' + {j^tdt -udu), 

3. a = n-4/("+2)^ n^2: 

^max ^ ^ker^^4^^^^(^^2)ua„, {2x'x^ - 6'^x\'^)dj - {n + 2)x'udu- {1 + lJ.){n + 2)x'vdy) , 

a = u-^, n = 2: A"^^^ = A^"' + {tdt + udu, (l){x\x'^)di +i>{x\x'^)d2 - 2(pi{x^ , x^)udu) , 
where (j){x^,x'^) and ipi^x^ ,x'^) are arbitrary solutions of the Cauchy-Riemann system (j)i = 

4. a = 1: 

^max ^ ^Q^^ g,^ ^ig, _ ^j q,^ a{t,x)dy, vdy, 2tdt + x'Oi, uOu, 2tdi - x'udu + x'vd,,, 
At^dt + Atx'di - (|xp + 2nt)udu + {\x\^ - 2nt)vdy, /3(t,x)9„), 
where a{t,x), j3{t,x) are arbitrary solutions of the linear heat equations at + an = d and 



It is the obvious extension of the known symmetry algebras [8, 9, 13] of diffusion equations to 
the adjoint equation. 

Applying the infinitesimal criterion for the variational invariance we get the basis of the 
algebra of variational symmetries 



A 



ker 



Pt = dt, Pi = di, Jij = x'dj - x^di, V{a{x)) = a{x)dy, Di = 2tdt + x'di - nvd^), 



of variational functional ()17() for arbitrary a{u) and its extension 

^iar = {Pt = dt, Pi = di, Jij = x'dj - x^di, Vi{a{t,x)) = a{t,x)dy, 

Di = 2tdt + x*9j — nvdv, D2 = udu — vd^, Gi = 2tdi — x^udu + x^vd^, 

n = 4t'^dt + Atx'di - (|xp + 2nt)udu + {\x\^ - 2nt)vd^, U{fl{t, x)) = (3{t, x)du, ) 

for a{u) = 1. For a{u) = e^, A^^, = A^% + {D2 = tdt - du). 

If a{u) = u^',^Ji^ 0, -4/(n + 2) then A^^, = A\Z + (^2 = ^^t^t - udu + vd^). 
For a = n-4/("+2), 



A. 



-- A^^'^ + {D2 = Atdt + (n + 2)udu, 
Ha = {2x^x^ - 5'^\x\^)dj - (n + 2)x'udu -{l + ^l){n + 2)x'vd^). 



According to the generalized Noether theorem there exists a one-to-one correspondence be- 
tween the found variational symmetries and the nontrivial conservation laws of system p6j) . 
Thus, using the well-known formulae for the components of the conserved vectors (see, e.g., [19, 
Corollary 4.30]) we obtain nontrivial conservation laws of system (|16|) presented in table 1. 

Table 1: Conserved vectors of nonlinear system H16|) . 



N 


Q 


T 


X' 


1 


V{a) 


a{x)u 


ctiA — aaui 






2 


Pt 


aUiVi 


—a{utVi + UiVt) 






3 


Pj 


\ (uVj — VUj ) 


5^^ {\(utv - uvt) + aujVj) - aUiVj 


~ aujVi 




4 


Jkl 


^v{uk ~Ul) + ^u{vi ~Vk) 


(uk - ui)avi + {vk - vi)aui 






5 


Di 


2taUjVj + ■^x-' (uVj — Ujv) + ^uv 


^x'{utv — uvt) — avi{x^Uj + 2tut) - 
+ax^UjVj 


- aui{nv + x-'vj 


+ 2tvt) + 



Here A — A{u) — J a{u)du, a — a{x) is an arbitrary solution of the Laplace equation an = 0, 5'-' = 1 if i = j 
and (5'-' = if i / j, 1 < i,j, k,l < n. 

System (|1()|) with a = e" admits extra conservation law of form Dt{te'^UiVi — v) — Di{t{utVi + 
UiVt) + avi) = 0. 

For case a = u^, /u 7^ 0, — 4/(n + 2), system (|l()j) has a CL with the conserved density and 
fluxes of form T = fitauiVi — uv, X'^ = —a{uvi — Uiv) — ^ta{utVi + UiVt). If /U = — 4/(n + 2) 
system (|lf)|) admits additionally 



Dt 
+D^ 



-{2x^x^ - 5''^\x\'^){ujV - UVj) - fj.{n + 2)x*t; 



(n + 2)(5*-^n"+2t) + (2x^x-' — 5^^x\ ) I -{utv — uvt) + a{ukVk — UjVi — UiVj) 



-{n + 2)x^a{uvi — (1 + fj)uiv) 



0. 



Conservation laws of the linear system are presented in table 2. 
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Table 2: Conserved vectors of system ut — uu = 0, vt + va = 



N 


Q 


r 


X' 


1 


Vr{a) 


au 


OiU — aui 


2 


Pt 


UiVi 


— {utVi +UiVt) 


3 


P, 


^{uVj — VUj) 


(5*^ (i(lltU — UVt) + UjVj) — UiVj — UjVi 


4 


Jki 


^v(uk - Ui) + ^u{vi - Vk) 


{Uk -Ul)vi + [vk -Vl)Ui 


5 


Di 


2taujVj + ^x^{uVj ~ Ujv) + ^uv 


\x^(utv — UVt) - Vi{x^Uj + 2tut) - Ui{nv + x^Vj + 2tvt) + 

+x'UjVj 


6 


[/(/?) 


I3v 


Pv, -(3,v 


7 


D2 


uv 


UVi ~ UiV 


8 


Gj 


— X-'UV — t{UjV + UVj) 


(5*^ [^{vut — UVt) + auiVi) — x^ (uVi + Uiv) — 2t{ujVi + UiVj) 


9 


n 


4f UiVi + 2tx^{uvi — Uiv) — \x\ uv 


2tx'{utV — UVt) + 4tx^UjVj — 4:tx^(ujVi + UiVj) — 
-uv^{\x\^ + 2nt) +Uiv{\x\^ -2nt) - 4i^{utVi + UiVt) 



Here a = a{t, x), P — j3{t, x) are arbitrary solutions of the linear heat equations ctt + an = and Pt — Pa, 5^-' = 1 
\i i = j and 5*-' = 1 if i 7^ j. 

Note that the conserved vectors from case 1 of tables 1-2 do not involve the solutions of the 
adjoint equation, and therefore, determine local conservation laws of diffusion equation (|15j) . 
These are exactly the conservation laws obtained above with application of the direct method 
(see theorem 121). The conserved vectors from cases 1-5 (1-9) of table 1 (table 2) determine the 
nonlocal conservation laws of diffusion equation (|15j) . where the non-local variable v is defined 
by the associated equation vt + avu = 0. 

Consider now in more detail the linear case 



Ut 



0, 



Vt + Vii = Q 



(19) 



Let us emphasize that the conserved vectors presented in table 2 are inequivalent with respect 
to the usual equivalence relation given in definition |31 and the corresponding conservation laws 
are linearly independent. At the same time, acting by the transformations from Lie symmetry 
group of system (|19|) . we can establish G^^r -equivalence relations on the space of conservation 
laws. Thus, case is equivalent to case ^ under the discrete symmetry transformation t = —t, 
X = X, u = —V, V = —u. 

Consider now the set of conservation laws spanned by those from cases ^ and [3 It is obvious 
that they form a linear subspace of the space of conservation laws of system ()19() with a set of 
generatrices {{au, aiu — aui), {uv, uvi — Uiv)} (as usual, we identify conservation laws with their 
representatives in the space of conserved vectors). Let us now change the set of generatrices 
as {{{v + au), {vi + ai)u — {v + a)ui), {uv,uvi — Uiv)}. Under the point transformation v — > 
V — a (the other variables are not transformed) from the symmetry group G\^-^ the conserved 
vector {{v + au), {vi + ai)u — {v + a)ui) is mapped to the {uv, uvi — Uiv). Therefore, up to the 
symmetry group of system (fT^ the considered subspace of conservation laws is generated by 
the conservation law with the conserved vector {uv,uvi — Uiv). 

In [12, 18] it is shown that for a system of Euler-Lagrange equations a generating set of in- 
equivalent conservation laws is determined by the ones with conserved vectors that correspond to 
symmetry operators spanning the symmetry algebra under the action of adjoint representations. 
Using this property we can complete investigation of equivalences between the conservation laws 
of system ^. Since [Pi,Gi] = D2, [Pt,Di] = 2Pt and [Pt,U] = ADi - 2nL>2, the conservation 
laws corresponding to Pt, Di and D2 are dependent on those generated by Pi, Gi, and 11. Using 
that [Pt, Gj] = Pj we deduce that the conservation law generated by Pj can be obtained from the 
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one corresponding to Gj under the multiple actions of symmetry transformations. It is obvious, 
that any conservation law from case |1] is equivalent to those generated by J12 with respect to the 
rotation transformations. Similarly, case|Slis G^3^j.-equivalent to the conservation law generated 
by Gi. Therefore, the following statement is true. 

Theorem 5. The generating set of conservation laws of system (|19|) consists of the conservation 
laws corresponding to the variational symmetry operators J12, Gi, H. 

Analogously we can find the generating set of conservation laws of nonlinear system (|l(i|) . 

Theorem 6. The generating set of conservation laws of nonlinear system (|19|) with arbitrary 
value of a{u) consists of the conservation laws corresponding to the variational symmetry oper- 
ators Di, J12 and V{a{x)). For o = e" and a = u^ all the conservation laws are generated by 
ones corresponding to Di, D2, J12 and V{a{x)). 

Note 2. There exist a non-trivial G^'^^-equivalence relation on the set of conservation laws of 
system ()19p corresponding to the variational symmetry operators V{a{x)), generated by the 
rotation transformations, translations of space variables and scale transformations. 

5 Conclusion 

In the presented work we used the direct method and the composite variational principle for 
construction of all possible inequivalent local conservation laws for class ^ of n-dimensional 
diffusion-convection equations. In contrast to the case of one-dimensional diffusion-convection 
equations, there exist multidimensional nonlinear equations with infinite-dimensional spaces of 
conservation laws. The composite variational principle allows us to determine a set of nonlocal 
conservation laws for the equations under consideration. 

The adduced results can be developed and generalized in a number of directions. Thus, 
e.g., the problem of the complete classification of conservation laws of the obtained system of 
Euler-Lagrange equations, which is closely related to construction of generalized (higher order) 
symmetry operators, is still open. Another promising direction is connected with investigation 
of potential symmetries and potential conservation laws. For such systems (see e.g. [1]) it follows 
from existence of gauge transformations for potentials that there exist no nontrivial potential 
symmetries without the introduction of a gauge. A natural question arises: How does one choose 
a gauge to obtain a nontrivial potential symmetry? (Existence of potential symmetries can be 
used as a criterion of adequacy of gauges applied in physical theories.) Moreover, we should like 
to classify gauges depending on group properties of potential symmetries up to symmetry group 
of the system under consideration. When potentials are introduced and a gauge is chosen, we 
can investigate symmetries of the potential system, construct its conservation laws and introduce 
new additional potential variables. For the extended potential system we again have a problem of 
choosing a gauge. Iterating the above procedure, we come to the problem of finding a hierarchy 
of right gauges. 
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